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ABSTRACT 



We present results of a high-resolution soft X-ray (0.2-2 keV) spectroscopic study 
of a sample of 69 nearby obscured Active Galactic Nuclei (AGN) observed with the 
Reflection Grating Spectrometer (RGS) on board XMM-Newton. This is the largest 
sample ever studied with this technique so far. The main conclusions of our study can 
be summarized as follows: a) narrow Radiative Recombination Continua are detected 
in about 36% of the objects in our sample (in 26% their intrinsic width is < 10 eV); b) 
higher order transitions are generally enhanced with respect to pure photoionization, 
indicating that resonant scattering plays an important role in the ionization/excitation 
balance. These results support the scenario, whereby the active nucleus is responsible 
for the X-ray "soft excess" almost ubiquitously observed in nearby obscured AGN via 
photoionization of circumnuclear gas. They confirm on a statistical basis the conclu- 
sions drawn from the detailed study of the brightest spectra in the sample. Further- 
more, we propose a criterion to statistically discriminate between AGN-photoionized 
sources and starburst galaxies, based on intensity of the forbidden component of the 
Ovil He-a triplet (once normalized to the Ovni Ly-a) coupled with the integrated 
luminosity in He- like and H-likc oxygen lines. 

Key words: galaxies: active - galaxiesmuclei - galaxies :Seyfert - X-rays: galaxies 



1 INTRODUCTION 

Nearby X-ray obscured Active Galactic Nuclei (AGN) in- 
variably exhibit excess emission above the extrapolation of 
the absorbed nuclear emission (Turner et al. 1997; Guainazzi 
et al. 2005). The origin of this component - which can rep- 
resent a significant fraction of the active nucleus bolometric 
energy budget (Levenson et al. 2002) - is still largely un- 
known. Gas in the nuclear environment could be heated to 
million degrees by shocks induced by AGN outflows (King 
2005) or episodes of intense star formation (Cid Fernandes 
et al. 1998; Gonzalez-Delgado et al. 2001). Alternatively, the 
AGN primary emission could photoionize and photoexcite 
circumnuclear gas. 

The latter scenario has recently received direct obser- 
vational support, thanks to high-resolution capabilities in 
the spatial and frequency domains that large X-ray observa- 
tories such as Chandra and XMM-Newton nowadays offer. 
High-resolution spectra unveiled signatures of photoionized 
plasma in a few bright objects: NGC 1068 (Kinkhabwala et 
al. 2002; Young et al. 2001), the Circinus Galaxy (Sambruna 
et al. 2001), Mkn 3 (Sako et al. 2000; Bianchi et al. 2005; 



Pounds et al. 2005); NGC 4151 (Schurch et al. 2004). In 
NGC 1068 the contribution of collisionally ionized plasma 
to the observed soft X-ray emission is constrained to be 
lower then 10% (Brinkman et al. 2002). These conclusions 
are based on three pieces of experimental evidence: 

(i) the spectra are dominated by strong emission lines of 
highly-ionized species from Carbon to Silicon, as well as by 
L-shell transitions from Fexvn to Fexxi 

(ii) narrow Radiative Recombination Continua (RRC) 
features from Carbon and Oxygen were detected. The width 
of these features indicates typical plasma temperatures of 
the order of a few eV (Kinkhabwala et al. 2002). These 
features are unequivocal signatures of photoionized spectra 
(Liedahl & Paerels 1996). 

(iii) the intensity of higher-order series emission lines, 
once normalized to the K Q , are larger than predicted by 
pure photoionization, and are consistent with an important 
contribution by photoexcitation (resonant scattering) (Band 
et al. 1990; Matt 1994; Krolik et al. 1995). This explains why 
standard plasma diagnostics (Porquet & Dubau 2000) fail to 
properly interpret the physical nature of the spectra. 
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A solution in terms of AGN-photoionized gas can also 
explain the coincidence in extension and overall morphology 
between soft X-ray emission and the Narrow Line Regions 
(NLRs), the latter traced by 0[lll] HST maps, on scales as 
large as a few hundred parsecs (Bianchi et al. 2006). Solu- 
tions satisfying the observed X-ray to optical flux ratio re- 
quire an approximately constant ionization parameter (i.e., 
a density scaling as the inverse square of the distance from 
the ionizing source), similarly to what is often found using 
photoionization models of the NLRs. 

So far, high-resolution X-ray spectra have been pub- 
lished only for a few X-ray bright obscured AGN. However, 
diagnostically important emission lines in these objects ex- 
hibit very large Equivalent Widths, EWs, as the contin- 
uum is often totally suppressed in the soft X-ray band. In 
this paper we present the first systematic high-resolution 
X-ray spectroscopic study on a sizable sample of obscured 
AGN. The Reflection Grating Spectrometer (RGS) on board 
XMM-Newton (der Herder et al. 2001) is the most suitable 
instrument currently flying for this purpose, due to its un- 
precedented effective area in the 0.2-2 keV band, as well as 
its good absolute aspect solution accuracy (~8 mA). 



2 CIELO-AGN 

2.1 The sample 

Our sample comprises all the type >1.5 AGN (according 
to the NED classification) observed by XMM-Newton, and 
whose data were public as of September 2006 (Tab. 1 shows 
the list of sources in the sample, together with their redshift 
and the combined exposure time for one RGS camera) . After 
excluding H1320+551 and H1419+480 due to uncertainties 
in their classification, the sample includes 69 sources. 

2.2 Data processing 

For each observation, we have reprocessed RGS data start- 
ing with the Observation Data Files, using SASv6.5 (Gabriel 
et al. 2003), and the most advanced calibration files avail- 
able as of May 2006. Background spectra were generated 
using blank field event lists, accumulated from different po- 
sitions on the sky vault along the mission. Spectra of the 
same source from different observations were merged, to- 
gether with their response matrices, after checking that no 
significant spectral variability occurred. 

Each spectrum was systematically searched for the pres- 
ence of emission lines. We simultaneously fit together the 
spectra of the two RGS cameras, using Xspec version 12.3.0. 
"Local" fits to the data were performed on the unbinned 
spectra on ~100 channels wide intervals, using Gaussian 
profiles to account for any line features. In only a few cases 
(see Tab. 2) the width of the profiles accounting for bound- 
bound transitions is inconsistent with zero; hence in these 
cases only the intrinsic profile width has been left free in the 
fit. On the other hand, free-bound transitions were modeled 
with Gaussian profiles, where the intrinsic width was always 
considered an additional free fit parameter. We have always 
assumed the same intrinsic width for the components of a 
multiplet. Local continua were modeled with F = 1 power- 
laws, leaving free in each fit the continuum normalization (r 



Table 1. CIELO-AGN sources. 



Source 


s 


Exposure time 






(ka) 


Circinus Galaxy 


0.001 


102.2 


ESO509-G66 


0.045 


13.6 


IC2560 


0.010 


81.4 


IC4395 


0.036 


22.1 


IC4995 


0.016 


11.5 


IIIZW03B 


0.027 


20.7 


IRAS01475-0740 


0.018 


11.5 


IRAS08572 + 3915 


0.058 


28.6 


IRAS09104+4109 


0.442 


13.5 


IRAS10214+4724 


2.286 


53.4 


IRAS13197-1627 


0.017 


44.6 


IRAS15480-0344 


0.030 


10.9 


MCG-5-23-16 


0.009 


61.6 


MRK3 


0.014 


180.8 


MRK6 


0.019 


96.4 


MRK331 


0.019 


17.9 


MRK348 


0.015 


46.3 


MRK612 


0.020 


11.8 


MRK744 


0.009 


22.7 


MRK993 


0.015 


22.6 


MRK1152 


0.053 


26.4 


NGC1068 


0.004 


86.8 


NGC1365 


0.005 


127.5 


NGC1386 


0.003 


16.9 


NGC1410 


0.025 


11.7 


NGC1614 


0.016 


23.5 


NGC2110 


0.008 


51.7 


NGC2273 


0.006 


12.5 


NGC2623 


0.019 


12.3 


NGC2992 


0.008 


28.5 


NGC34 


0.020 


21.8 


NGC3982 


0.004 


34.0 


NGC4138 


0.003 


14.3 


NGC4151 


0.003 


55.5 


NGC4168 


0.007 


23.0 


NGC424 


0.012 


7.9 


NGC4258 


0.002 


78.5 


NGC4303 


0.005 


42.8 


NGC4395 


0.001 


108.4 


NGC4472 


0.003 


1.5 


NGC4477 


0.004 


13.3 


NGC449 


0.016 


11.9 


NGC4507 


0.012 


44.5 


NGC4565 


0.004 


14.4 


NGC4639 


0.003 


14.6 


NGC4725 


0.004 


17.6 


NGC4945 


0.002 


63.8 


NGC4968 


0.010 


11.8 


NGC5033 


0.003 


19.7 


NGC5252 


0.023 


65.8 


NGC526A 


0.019 


46.9 


NGC5273 


0.004 


16.1 


NGC5506 


0.006 


32.9 


NGC5643 


0.004 


9.6 


NGC591 


0.015 


11.8 


NGC6552 


0.026 


6.8 


NGC7172 


0.009 


72.0 


NGC7212 


0.027 


13.7 


NGC7314 


0.005 


43.4 


NGC7479 


0.008 


12.6 


NGC7582 


0.005 


22.3 


NGC7674 


0.029 


10.3 


UGC1214 


0.017 


11.8 


UGC2456 


0.012 


17.1 


UGC2608 


0.023 


7.5 


UGC4203 


0.014 


7.8 


UGC6527 


0.027 


23.6 


UGC8621 


0.020 


11.8 


UM625 




11.6 



is the power-law photon index) . No assumption was made a 
priori on the line centroid energies. 1 However, fits of He-a 
triplets have been performed keeping the relative distance 
between the centroid energies of the components fixed to 
the value dictated by atomic physics. A line is considered to 
be detected when its flux is inconsistent with at the l-cr 
level. In Appendix A we list energies, fluxes and intrinsic 
widths (for the RRC only) of the transitions individually 
presented in this paper, and discuss the consistency of our 
measurements with the expected laboratory energies [labo- 
ratory energies are extracted from the CHIANTI database 
(Dere et al. 2001)]. Line luminosities have been corrected 
for Galactic photoelectric absorption using column densities 
after Dickey & Lockman (1990). 



1 On the other hand, upper limits on the line intensities were 
calculated assuming the laboratory energies. 



On the origin of soft X-rays in obscured AGN 3 



Table 2. 1-ct intrinsic widths for resolved Gaussian profiles in 
CIELO-AGN. 



Transition 


Width 




(km s- 1 ) 


NGC1068 



Fcxvn 3d-2p (iPi) 
Nvn Ly-a 
Ovn He-a 
Ovu He-/3 
Ovn He-7 
Ovm Ly-a 
Ovm Ly-/3 



670±JJ° 
820±f° 
550±5° 
3300±1»» 
650±i5° 
460±|0 

l_200 



1000± 



20 



NGC4151 
Fexvm 3s-2p 
Cvi Ly-a 
Cvi Ly-/3 



llon+ 900 

490± 80° 
nnn+700 



NGC4507 
Ovn He-a 



520±i 



We have as well checked whether the detected line en- 
ergy centroids in each source were systematically shifted 
with respect to the laboratory energies. We define "system- 
atic" energy shifts different from zero (in the same direction) 
at the lcr level in at least two of the following transitions: 
Cvi Ly-a, Ovn He-a, Ovn He-/3, and Ovm Ly-a. Our anal- 
ysis found none of these systematic shifts. Nonetheless, we 
found some occurrences of discrepancies between the best-fit 
centroid and the laboratory energies in individual lines. We 
refrain from attributing any astrophysical meaning to these 
discrepancies, which may be due to contamination by nearby 
transitions, residual errors in the aspect solution of the RGS 
cameras, or spurious detections. The reader may find a more 
detailed discussion of this point in the Appendix. 

In Fig. 1 we show the RGS spectra of the three still 
unpublished Seyfert 2 galaxies in our sample, which exhibit 
the largest number of line detections. 

Hereafter, uncertainties on the fitting parameters are 
quoted at the lcr level; likewise the upper limits represent 
a lcr confidence level. Errors on the line centroid energies 
include a 8 mAsystematic error. Throughout this paper, en- 
ergies/wavelengths of astrophysical lines are quoted in the 
source rest frame, unless otherwise specified. 



2.3 The catalog 

The results of the procedure outlined above have been 
compiled into CIELO-AGN ( Catalog of Ionized Emission 
Lines in Obscured AGN). Our study exhibit a high-detection 
efficiency, at least with respect to the most common lines ob- 
served in obscured AGN spectra, despite the overall low soft 
X-ray flux (sample median 6xl0 -13 erg cm~ 2 s _1 ; see Tab. 3 
for a complete list of the lines detected in CIELO-AGN). 



3 RESULTS 

In this Section we will use some results of our study to try 
and answer the following questions on the nature of soft 
X-ray emission in obscured AGN: 



Table 3. Number of CIELO-AGN sources, N det , for which a 
given Transition (laboratory wavelength A; a (,) has been detected, 
if N det > 



Transition 


(A) 


Ndet 


Sixm He-a 


6.740 


5 


Mgxi He-a (r) 


9.228 


12 


Mgxi Hc-a (/) 


9.314 


13 


Nex Ly-a 


12.134 a 


20 


Fcxxi 3d-2p 


12.282 


21 


Fexx 3d-2p 


12.845 


19 


Neix He-a (r) 


13.447 


13 


Neix He-a (i) 


13.553 


6 


Neix Hc-a (/) 


13.699 


15 


Ovm RRC 


14.228 


13 


Fcxvni 3d-2p 


14.413'' 


24 


Ovm Ly-8 


14.821 c 


20 


Fexvn 3d-2p ( x Pi) 


15.015 


25 


Fexvn 3d-2p ( 3 £>i) 


15.262 


24 




1D.UUO 


99 


Fexvm 3s-2p 


16.091 


23 


Ovn RRC 


16.771 


17 


Fcxvn 3s- 2p (3G/M2) 


17.076 d 


27 


Ovn He-7 


17.768 e 


16 


Ovn He-/3 


18.627 


20 


Ovm Ly-a 


18.967 


30 


Ovn Hc-a (r) 


21.602 


22 


Ovn He-a (i) 


21.803 


13 


Ovn He-a (/) 


22.101 


23 


Ncvn Ly-a 


24.782^ 


24 


Cvi Ly-/3 


28.459 s 


14 


Nvi He-a 


28.787 


11 


Cv RRC 


31.622 


6 


Cvi Ly-a 




23 


Ai = 12.1321 A, A 2 = 12.1375 A 



^triplet: Ai = 14.3760 A, A 2 = 14.4187 A, A 3 = 
c possible contamination by Fexx 3p 2 -2p 2 at A 



14.4210 A 
= 14.8501 A 



at A = 17.8472 A 



d doublet: Ai = 17.0500 A, A 2 = 17.0970 A 
e possible contamination by Fexvm 3p 2 — 2p 2 
f doublet: \± = 24.7790 A, A 2 = 24.7840 A 
^doublet: Ai = 28.4612 A, A 2 = 28.4662 A 
^doublet: Ai = 33.7342 A, A 2 = 33.7396 A 



• is AGN photoionization the dominant physical process? 

• does resonant scattering play an important role? 

• do efficient X-ray line diagnostics exist, which could al- 
low us to discriminate on a statistical basis between AGN- 
and starburst-powered spectra? 

3.1 RRC 

The three most intense RRC transitions in the RGS energy 
bandpass are: Ovm at 14.228A, Ovn at 16.77lA and Cv at 
31.622A. At least one of these features is detected in 36% 
(25/69) of the objects of our sample. The few measurements 
of the RRC width, which represents a direct estimate of the 
temperature of the plasma (Liedahl & Paerels 1996), cluster 
in the range 1-10 eV. There is an obvious selection effect, 
favoring the detection of narrow features. Nonetheless, still 
in 29% (26%) of the objects in our sample the RRC width 
is constrained to be lower than 50 (10) eV (see Fig. 2). In 
10 further objects we detect upper limits on the RRC lumi- 
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Figure 1. RGS spectra for the three galaxies of our sample, which exhibit the largest number of detected emission lines (in brackets 
below) and whose RGS spectra have not been published yet: NGC 1365 (19), NGC 4507 (17), UGC 1214 (14) [for comparison the number 
of detected lines in the brightest Seyfert 2 of our samples are: 25 (NGC 1068), 19 (the Circinus Galaxy), 18 (NGC 4151), 17 (Mrk 3)]. 
Spectra of the two RGS cameras have been merged and smoothed with a 5-channcls wide triangular kernel for illustration purposes only. 
The positions of the line transitions measured in CIELO-AGN are labeled. 
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Figure 2. Luminosity versus intrinsic width for the RRC fea- 
tures detected in our sample. Filled data points correspond to a 
measurement of the RRC width; empty data points correspond to 
RRC width upper limits. Circles: Ovill; squares: Ovn; triangles: 
Cv. Data points corresponding to upper limits on both quantities 
are not shown for clarity. 
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CVII RRC (16.777A)/FeXVll 3s-2p (17.076A) intensity ratio 



A 



1 2 3 4 5 6 7 

Figure 3. Observed intensity ratio between the 16.777A and the 
17.076A features in CIELO-AGN. Shaded cells indicate lower lim- 
its. The dashed line represents the maximum expected ratio, if 
both features are due to the Fexvn triplet. In 15 out of 16 cases, 
the observed ratio exceeds this limit, indicating that Ovn must 
be the dominant transition at 16.777A. 



nosity in the range where measurements are found: in 7 (10) 
of them the luminosity is constrained to be < 10 41 ' 1 erg s _1 
(< 10 44 ' 2 erg s _1 ). Finally, in 34/69 (49%) of the objects 
the quality of the data is too poor to allow a significant de- 
tection of any of the RRC features, and the upper limits on 
their luminosities are inconclusive. 

In principle, measurements of the Ovn RRC could be 
contaminated by the 3F component of the Fexvn triplet at 
16.780A (B rown et al. 1998). However, the intensity ratio 
between the Ovn RRC and the 3G/M2 components of the 
same triplet (which cannot be resolved by the RGS) exceeds 
0.6, the largest expected value if the former feature were en- 
tirely due to the 3F component (Phillips et al. 1997; Mauche 
et al. 2001; Beiersdorfer et al. 2002), in 15 out of 16 cases 
(by more than a factor of 2 in 13 out of 16; cf Fig. 3). We 
are therefore confident that most of our measurements of the 
the emission line at 16.777A represent bona fide detections 
of the Ovn RRC. The total number of sources for which at 
least one RRC feature is detected is reduced by 1 if the effect 
of this potential mis-identification is taken into account. 

3.2 Higher order series 

In addition to n = 2p-^ls transitions for H- and He-like 
atoms, CIELO-AGN contains a fair number of detections of 
discrete higher order resonance transitions (np— >ls, n > 2). 
These transitions are selectively enhanced by photoexcita- 
tion. Since the forbidden / transition in He-like triplets is 
unaffected by photoexcitation, the intensity ratio between 
higher order series and / transition intensities provides a 
potentially powerful diagnostic of the importance of reso- 
nant scattering in radiation ionized spectra (Kinkhabwala 



et al. 2002). An example of the application of this diagnos- 
tic test to CIELO-AGN is shown in Fig. 4, where we display 
the intensity of the Ovn He-/3 against the / component of 
the Ovn He-a triplet. We compare the experimental results 
with the predictions of pure photoionization, and with mod- 
els where radiative decay from photoexcitation and recombi- 
nation from photoionization are self-consistently calculated 
(model photoion; Kinkhabwala et al. 2002). We have pro- 
duced a grid of models for different values of Ovn column 
densities (Novn G[10 15 ,10 20 cm -2 ]), turbulence velocities 
{vturb €[0,500 km s -1 ]), and temperatures (kT G[l,20 keV]). 
The weighted mean of the He-/3 versus / intensity ratio is 
0.25 ±0.03, larger then expected for pure photoionization. It 
corresponds to Ovn column densities Novn — 10 17-18 cm -2 
(the dependence on the other parameters is small). The mea- 
surements on the He-/3 Ovn transition at 18.6270A could be 
contaminated by the nearby Nvn RRC at 18.5872A 2 . In or- 
der to minimize this contamination, we fit the X-ray spectra 
around the Ovn He-/3 feature in a range, which in principle 
does not include the contaminating feature. Moreover, we 
have verified that similar enhancements of higher-order lines 
to the Lyman-a of H-like Carbon and Oxygen are observed 
as well (cf. Tab. 4). The dependence on the temperature, col- 
umn density and velocity distributions in H-like species is 
such that no meaningful constraints on column density can 
be drawn. However, in 8 of the individual brightest sources 



2 For instance, the Ovn He-/3 width measured in in NGC 1068, 
a = 73±| 3 eV (cf. Tab. 2), the largest measured in this spectrum, 
is most likely contaminated by the Nvn RRC 
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Ratio between OVII He-/? and forbidden (f) 




10" 7 10" 6 10 s 10" 4 

OVII He-/S Intensity (ph cm 2 s 1 , 18.616A) 







(ks) 


ARP244 


0.006 


75.0 


ARP270 


0.005 


36.6 


Holmbergll 


0.000 


30.6 


IZwlS 


0.003 


52.3 


M82 


0.001 


141.6 


M83 


0.002 


30.5 


MOG- 5-31-7 


0.010 


19.0 


NGC1313 


0.002 


241.4 


NGC1482 


0.006 


17.9 


NGC1511 


0.004 


44.0 


NGC1705 


0.002 


58.5 


NGC2146 


0.003 


35.9 


NGC2403 


0.000 


43.5 


NGC253 


0.001 


192.3 


NGC2798 


0.006 


13.8 


NGC3256 


0.009 


16.7 


NGC3310 


0.003 


19.2 


NGC3690 


0.010 


21.5 


NGC4214 


0.001 


19.5 


NGC4449 


0.004 


28.1 


NGC474 


0.008 


25.4 


NGC5073 


0.009 


51.8 


NGC520 


0.008 


12.6 


NGC5253 


0.001 


45.8 


NGC660 


0.003 


11.9 


NGC7552 


0.005 


24.0 


NGC7714 


0.009 


41.3 



Table 5. Log of the XMM-Newton observations of the starburst 
galaxies in the control sample 



Figure 4. Intensity of Ovil He-/3 line against the intensity of the 
/ component of the He-a triplet (only data points correspond- 
ing to a detection of the latter are shown; data points corre- 
spond to upper limits on the intensity of the former are shown 
as empty symbols). The dashed- dotted lines represent the predic- 
tion of the photoion code for Ovn column densities increasing 
from Nqvii = 10 15 to 10 20 cm -2 in steps of one decade, assum- 
ing kT = 5 eV and v tur i, = 200 km s — 1 . The long-dashed line 
represents the predictions for pure photoionization. The shaded 
areas represent the loeii of the photoion predictions, when the 
turbulent velocity varies in the range 0-500 km s _1 at constant 
temperature and for the extreme values of the column density 
interval. The areas representing a variation of the temperature in 
the range 1-20 keV at constant velocity and column density are 
comparatively smaller, and therefore not shown. 



Table 4. Higher-order intensity ratios for selected transitions in 
CIELO-AGN. Each ratio is calculated as the weighted mean of the 
individual ratios on the N sources, where a measurement of both 
transitions is available, and using the statistical uncertainties on 
the ratio as weights. Expected values for pure photoionization 
(PIE) and collisional ionization (CIE) arc extracted from Tab. 3 
in Kinkhabwala et al. (2002) 



Ion 


Ratio 


TV 


PIE 


CIE 


Cvi Ly-/3/Ly-a 


0.51 ±0.12 


9 


0.14 


0.09 


Ovn He-/3// 


0.25 ±0.03 


12 


0.05 




Ovn He-j/f 


0.20±0.09 a 


9 


0.017 




O vni Ly-/3/Ly-a 


0.35 ±0.05 


1.3 


0.14 


0.10 



a this measurement is likely to be affected by contamination of 
Fcxviii 3p 2 -2p 2 at A = 17.8472 A 

in our sample 3 the ratio between the Ovm Ly-/3 and Ly-a 
intensities is large enough to be formally inconsistent with 



collisional ionization, and in 4 of them 4 with pure photoion- 
ization as well. 

3.3 AGN and starburst soft X-ray spectra 

We have analyzed the RGS spectra of a sample of 27 Star- 
Burst (SB) galaxies extracted from the Wu et al. (2002) 
sample (cf. Tab. 5) and observed by XMM-Newton, to iden- 
tify diagnostic criteria, which may allow us to statistically 
discriminate between SB- and AGN-dominated soft X-ray 
spectra. This sample has been analyzed in the same way 
as the Seyfert 2 sample. Some of the criteria quoted in the 
literature provide inconclusive results. For instance, the dis- 
tribution functions of the luminosity ratio between L shell 
iron and integrated oxygen He- and H-like lines are indistin- 
guishable (Fig. 5). Kallman et al. (1996) had already pointed 
out that L-shell iron transitions can give an important con- 
tribution to the overall luminosity budget in photoionized 
nebulae. The standard G ratio (G = (/ + i)/r; Gabriel & 
Jordan 1969; Porquet et al. 2000), is also ambiguous. Previ- 
ous studies have shown (cf. Sect. 1) that the enhancement 
of resonant lines by resonant scattering mimics the behavior 
of a collisionally ionized plasma. In the collisionally ionized 
case, the fluxes of all the triplet lines are enhanced, with the 
greatest increase in the resonance line. Resonant scattering 
boosts all of the resonant lines with respect to the forbidden 
and intercombination lines in the triplet. The value of G is 
therefore decreased by both collisional excitation and res- 
onant scattering. Finally, the fraction of SB spectra where 
narrow RRC features are detected is comparable to that 
observed in the Seyfert 2 sample, indicating that photoion- 
ization plays an important role in the ionization balance of 
starburst galaxies as well. 

In Fig. 6 we compare the distribution of the Ovn / 
transition intensity, normalized against the intensity of the 
Ovm Ly-a, as a function of the total luminosity in oxy- 
gen lines (integrated on all He- and H-like transitions). 



3 the Circinus Galaxy, IRAS 13197-1627, Mrk 3, NGC1068, 
NGC1365, NGC2110, NGC2992, NGC4945 



4 NGC1068, NGC1365, NGC2110, NGC2992 
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log (L Fe /L ) 
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Figure 5. Distribution function for the logarithm of the ratio 
of the integrated luminosities in the iron and oxygen lines for 
our Seyfert 2 (lower panel) and in the control sample of star- 
burst galaxies (lower panel). Diagonally shaded cells indicate up- 
per limits, checked cells indicate lower limits. The medians of the 
two distributions are very similar: —0.10 and —0.08, respectively. 



Ratio between OVIII K a and OVII K a forbidden (f) 
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Figure 6. Intensity of the / component of the Ovn triplet (nor- 
malized to the Ovm Ly-a intensity) against the total luminos- 
ity in Oxygen lines. Empty circles represents the obscured AGN 
in GIELO-AGN, filled circles the control sample of 27 starburst 
galaxies. Symbols representing censored data are plotted with a 
smaller size for the sake of clarity. 



luminosity is ~10 erg s~ in the former and ~10 erg s~ 
in the latter (using strict detections only). Still, the overlap 
between the line diagnostic distributions is significant, and 
prevents strong statements on individual sources. This most 
likely reflects the intrinsic "composite" nature of several ob- 
scured AGN (Cid Fernandes et al. 2001). 



4 CONCLUSIONS 

In this paper we present results from a systematic study of 
high-resolution soft X-ray spectra for 69 obscured AGN ob- 
served with the XMM-Newton RGS. Our main conclusions 
can be summarized as follows: 

- Radiative Recombination Continua: we detect RRC 
features in about 40% of the sample sources, and in 33% 
(27%) they allow us to constraint the temperature of the 
plasma < 50 eV (< 10) eV. This indicates that in a large 
fraction of objects in our sample photoionization dominates 
the ionization balance. Still, these percentages represent 
probably lower limits to the true fraction of photoionized 
sources in our sample: upper limits on narrow RRC features 
in almost 50% of the objects in our sample are still incon- 
clusive due to the lack of statistics 

- Resonant scattering: higher order transitions are en- 
hanced with respect to the expectation of pure photoion- 
ization, and inconsistent with the expectation of collisional 
ionization as well. This indicates that resonant scattering 
plays an important role in the ionization/excitation bal- 
ance. The observed Ovn He-/3 versus / intensity ratios 
are consistent with Ovn column densities in the range 
Novu ~10 17 ~ 18 cm" 2 . Interestingly enough, these values 
are in good agreement with the column densities measured 
in "warm absorbers" observed along the lines of sight to un- 
obscured AGN (Blustin et al. 2003; Steenbrugge et al. 2005; 
Blustin et al. 2005). 

- Starburst contribution: the comparison between the 
spectra of obscured AGN in our sample and a control sam- 
ple of nearby starburst galaxies suggests two empirical crite- 
ria to discriminate on a statistical basis between AGN- and 
starburst-powered spectra: a) total Oxygen line X-ray lumi- 
nosity <; 10 40 erg s _1 ; b) ratio between the Ovn triplet / 
component and the Ovm Ly-a <;1. 

The currently available X-ray instrumentation allowed 
us to explore the nature of the soft X-ray emission in ob- 
scured AGN as weak as ~0.03 mCrab. From this study, we 
have gained confidence that conclusions on the properties of 
the gas in the circumnuclear region of AGN extracted from 
the detailed study of high-quality spectra of the brightest 
objects can be extended to the whole population of nearby 
obscured AGN. This is the main message we would like to 
convey with this paper, which opens interesting perspectives 
for future enlargements of this study once deeper exposures 
or more sensitive instrumentation are available. 



AGN are generally characterized by larger intensity ratios, 
V- Vagn = 1.38±g;l|, r]SB = 0.57±8:o9 (5.8cr difference; the 
upper limits have been taken into account with a "boot- 
strap" method for censored data as in Schmitt 1985). The 
same figure also shows that the average line luminosity in 
AGN is larger than in starbursts. The median Oxygen lines 



APPENDIX 

The tables in this Section list centroid energies and fluxes for 
the bound-bound transitions individually discussed in this 
paper: Cvi Ly-a, Cvi Ly-/3, Ovn He-a, Ovn He-/3, Ovm 
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Source 


CVI Ly-a 


CVI Ly-/3 


OVII Hc-a (r) 


Ovu He-/3 


OviII Ly-a 


OviII Ly-/3 




(A, 33.737) 


(A, 28.459) 


(A, 21.602) 


(A, 18.627) 


(A, 18.967) 


(A, 16.006) 



Circinus Galaxy 

ESO509-G66 

IC2560 

IC4395 

IC4995 

IIIZW035 

IRAS01475-0740 

IRAS08572 + 3915 

IRAS09104+4109 

IRAS10214+4724 

IRAS13197-1627 

IRAS15480-0344 

MCG-5-23-16 

MRK3 

MRK6 

MRK331 

MRK348 

MRK612 

MRK744 

MR.K993 

MRK1152 

NGC1068 

NGC1365 

NGC1386 

NGC1410 

NGC1614 

NGC2110 

NGC2273 

NGC2623 

NGC2992 

NGC34 

NGC3982 

NGC4138 

NGC4151 

NGC4168 

NGC424 

NGC4258 

NGC4303 

NGC4395 

NGC4472 

NGC4477 

NGC449 

NGC4507 

NGC4565 

NGC4639 

NGC4725 

NGC4945 

NGC4968 

NGC5033 

NGC5252 

NGC526A 

NGC5273 

NGC5506 

NGC5643 

NGC591 

NGC6552 

NGC7172 

NGC7212 

NGC7314 

NGC7479 

NGC7582 

NGC7674 

UGC1214 

UGC2456 

UGC2608 

UGC4203 

UGC6527 

UGC8621 

UM625 



33.430±°;?JS> 



33 415 + - 074 
33.4±O=t 



33.767± 
33.776± 



0.052 

m 

0.064 



33.733± 
33.617± 



0.118 

:??7 

0.148 



33.895±°;°«« 



33 819 + 039 
0.060 
3d m 8 + 0. 123 



33.676±°°g 

33.545 + 0-057 
33.595±B;»SS 

33.797±°;J?| 



28.400±2'?2| 

28 43Q + o:095 
2O.439± 04g 



28.506 + 0-047 



28.152 + 0-086 



21.689 + 0-^4 
0.560 



28.491±[J[j|| 21.600±[ 



029 
522 



21.563 + 0-658 



91 61 9 + " -552 
21.012± 



21.602 + " -J 6 ! 



21.601 + 



0.029 



21.609 + 
21.573± 



0.530 
0.058 



1 0.049 
c 0.048 



18.367+0-392 



18.685 + 
18.634 + 



0.044 

m 

0.054 



18 740 + 092 



18.651 + ""- 
18.719±»» 4 8 



,78:91? 

"0.082 



1 0.040 
E 0.046 



18 931 +0.050 
18.931± 050 



19.013 + 0-J15 



18.962±[ 
18.966± 



lit: 

053 



18.813 + 001 
19.245±0»^ 



049 



18.992 + 0-072 



18.974±°- D3S 
18.986 + 8:813 



19.025 + 0-048 

■m 

0.118 



18.896 + 
18.928± 



18 976 + 0- 053 
0.052 



18.974 + 
18.977 + 



18.868 + 0-080 
18.948 + 0-OSS 

is Q«n+o:oS3 



ie.oiB±g : g*g 



1 6 099 + .061 
16.U22± 04g 



Mi 



15.961±0-;87 

15 - 819 ±o:q1q 

16.078 + 0-0S3 



I'M 



16.074± 
16.178 + 



16 001 + D11 ° 



16.065 + 0-141 
15.836±» : 4 44 

16.022 + 0-J03 



Table 6. Best-fit wavelengths for the transitions discussed in this paper. Dots indicates either lines not detected, or lines for which no 
constraints can be obtained (the luminosities in Tab. 7 are calculated in this case assuming the laboratory energies). In the table header, 
the laboratory wavelengths are in brackets. 



Ly-a, and Ovm Ly-/3 (Tab. 6 to 7). The same quantities 
are listed for the Cv, Ovu, and Ovm RRCs, alongside their 
intrinsic width (Tab. 8 to 9). 

The comparison between the best-fit centroid energy of 
detected lines in CIELO-AGN and the laboratory energies is 
an instructive exercise on the reliability of the catalog mea- 
surements. At the 2a level, the number of He- and H-like 



Oxygen lines (out the total detected ones) which are incon- 
sistent with the laboratory energies are: 0/25, 2/20, 2/30 and 
5/22 for Ovu Hc-a, Ovu Hc-/3, Ovm Ly-a and Ovm Ly-/3, 
respectively. The shifted measurements in the last case clus- 
ter around 783 eV (4 out of 5), where a strong Fexvin 3s2p 
transition is located (A = 15.83 A), and are therefore proba- 
bly due to a mis-identification. If we remove these suspicious 
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Cvi L 
10- 5 ph ci 



Circinus Galaxy 

ESO509-G66 

IC2560 

IC4395 

IC4995 

IIIZW035 

IRAS01475-0740 

IRAS08572 + 3915 

IRAS09104+4109 

IRAS10214+4724 

IRAS13197-1627 

IRAS15480-0344 

MCG-5-23-16 

MRK3 

MRK6 

MRK331 

MRK348 

MRK612 

MRK744 

MRK993 

MRK1152 

NGC106S 

NGC1365 

NGC1386 

NGC1410 

NGC1614 

NGC2110 

NGC2273 

NGC2623 

NGC2992 

NGC34 

NGC3982 

NGC4138 

NGC4151 

NGC4168 

NGC424 

NGC4258 

NGC4303 

NGC4395 

NGC4472 

NGC4477 

NGC449 

NGC4507 

NGC4565 

NGC4639 

NGC4725 

NGC4945 

NGC4968 

NGC5033 

NGC5252 

NGC526A 

NGC5273 

NGC5506 

NGC5643 

NGC591 

NGC6552 

NGC7172 

NGC7212 

NGC7314 

NGC7479 

NGC7582 

NGC7674 

UGC1214 

UGC2456 

UGC2608 

UGC4203 

UGC6527 

UGC8621 

UM625 



98.6±f<! 
1 1±0 3' 
H 



41 6±^ 9 
32.8 
0.6±H 



9.3±i 2 n 5 



ViA) 



1 13.9 
c 9.6 



10.9±}j> 
4±<" 



0.8 
0.7 



2.0±°f 



n 94-1-2 



,1.4 
c 1.3 



0.8±0| 
0.7± 



2.2 
1.1 



1 48.2 
c 32.0 



124.9±|-'5 
1.4±t>S 

18 ±o:l 



1 2+ 1 - 5 
1 4± 1 -3 

49.5±|;3 

05± R:a 

S.OiJi 
8.0±2-S 

0.9±!» 

i.9±«:| 
2.o±?;° 



2 2± b 



°- 3± o 

0.8±9 



1.0± 



1.0 
0.6 



,3.4 
c 2.5 



4±°' 6 



44.3±j| 



0.4±0;5 



1 3 + °' 8 
n 0.2 



0.2±° 



1 9 + °' 6 

o.2±°J 



14 9 + 7 - 3 



0.5 
0.4 



53.2±^.l 



1.9± 
1.0± 



M 
0.5 



15.5±1;3 



0.6± 
0.7± 



0.8±0;6 

°- 9 ±o:l 



20± xi 

2.0±2;3 
0.3±0-§ 



2.6±J; = 



0.7± 
2.3 



?:8 

+8:9" 
^0.7 



t 2.1 
c 1.5 



1 0.3 
c 0.2 



,0.2 
c 0.2 



14.9±J ■? 



,1.0 

c 0.8 



2 5±°' 6 
0.3±0 ; 3 



,8:1 



,¥:? 



0.4±0.3 



0.8± 
1.1± 



1.2 

8:0 
0.7 



,2.0 
c 1.5 



Table 7. Best-fit fluxes for the transitions discussed in this paper. Dots indicate lines not detected. 



cases, the ratio between the He-/? and the / transition in- 
tensity in the remaining objects is still inconsistent with the 
expectation of pure photoionization: 0.26 ±0.04 (cf. Tab. 4). 
The situation for the Cvi lines is more controversial. The 
centroid energies of 6 out of 23 Ly-a and 5 out of 14 Ly-/3 
measurements are inconsistent at the 2a level with the lab- 
oratory energies. There is no obvious explanation for these 
systematic discrepancies. No other potentially contaminat- 
ing transition is present, which could confuse the identifi- 



cation of the Cvi lines. Moreover, the shifted measurements 
do not cluster in any well defined spectral region. Uncertain- 
ties on the aspect solutions can be as high as 30A in this 
wavelength range, still insufficient to explain the observed 
differences. 

A similar analysis on the centroid energies of the RRC 
features shows that only 3 measurements (out of 16) of the 
Ovu RRC are inconsistent with the laboratory energy (at 
the 2a level; none for the Cv or Ovm RRCs). In this case, 
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Source 


Cv RRC 




Ovn RRC 


Ovm RRC 




E c 




E c a 


E c a- 




(A, 31.622) 


(eV) 


(A, 16.771) (cV) 


(A, 14.228) (cV) 



Circinus Galaxy 

ESO509-G66 

IC2560 

IC4395 

IC4995 

IIIZW035 

IRAS01475-0740 

IRAS08572 + 3915 

IRAS091044-4109 

IRAS10214+4724 

IRAS13197-1627 

IRAS15480-0344 

MCG- 5-23-16 

MR.K3 

MR.K6 

MRK331 

MRK348 

MRK612 

MRK744 

MRK993 

MRK1152 

NGC1068 

NGC1365 

NGC1386 

NGC1410 

NGC1614 

NGC2110 

NGC2273 

NGC2623 

NGC2992 

NGC34 

NGC3982 

NGC4138 

NGC4151 

NGC4168 

NGC424 

NGC4258 

NGC4303 

NGC4395 

NGC4472 

NGC4477 

NGC449 

NGC4507 

NGC4565 

NGC4639 

NGC4725 

NGC4945 

NGC4968 

NGC5033 

NGC5252 

NGC526A 

NGC5273 

NGC5506 

NGC5643 

NGC591 

NGC6552 

NGC7172 

NGC7212 

NGC7314 

NGC7479 

NGC7582 

NGC7674 

UGC1214 

UGC2456 

UGC2608 

UGC4203 

UGC6527 

UGC8621 

UM625 



31 486 + - 034 1 6 + °- 2 

31.511± 1 - b± 0.1 



31.626±0 : ;|° 



16 747+O O 42 
±0. ' 4 '=l=o.o51 



16.817± 
16.775± 
16.778± 



16.775± 
16.730± 



0.056 
I 

0.041 



16.929±°-°jjj- 
1 7 033 + 0. 138 



17 041± 0126 



< 33.2 

< 10.0 



9 o 1 y.B 

2ci± 1.0 



< 2.2 
2 1 + - 7 

< 13.9 



< 5.6 
< 22.0 



2 2+ L9 14 166 + 069 



14 197 + l ''°^ ;i 
" .044 



14.197±J 
14.453± 
14.234±; 



061 

m 

oso 



1 0.061 
E 0.051 



< 100.7 

2 - 7 ±o:8 
< 1.9 



12.4±8; 2 

< 66.6 
7-0-iV 



i.8±3 : ; 

< 5.2 



Table 8. Best-fit wavelengths and widths for the RRCs discussed in this paper. Dots indicates either lines not detected, or lines for which 
no constraints can be obtained (the luminosities in Tab. 9 are calculated in this case assuming the laboratory energies.) 



as discussed in Sect. 3.1, the discrepancy is probably due to 
confusion with the nearby Fexvn triplet. 
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Cv RRC Ovn RRC Ovm RRC 

If)- 5 ph cm- 2 s- 1 lO" 5 ph cm" 2 s- 1 pi, en,- 2 



Circinus Galaxy 




~ . 1 1 9 3 
2 4 ±2 3 


„ „ 1 1 7 


ESO509-G66 








IC2560 




0.3±0| 




IC4395 






IC4995 








IIIZW035 








IRAS01475-0740 








IRAS08572 + 3915 








IRAS 09 104+4109 








IRAS10214-j-4724 






0<i± 03 


IRAS13197-1627 






IRAS15480-0344 






MCG-5-23-16 




•** 




MRK3 




1 6± o.l 


" 5 

°- 9± o:3 


MR.K6 








MRK331 








MRK348 








MRK612 








MRK744 








MRK993 








MRK1152 






., „ „ . 27 4 
10.1±f ' 


NGC1068 


„ , in 1 


65. 9± 2 '} 


MPfi QfiC 
1\Kj<^j LOOO 


67. 2±i ^ 

°- 8± c>:5 


1 i+W 


n .s + 0.2 


NGC1386 


°- 3 ±o:!j 


NGC1410 






NGC 1 6 14 








NGC2110 






6 

1,4± o!5 


NGC2273 






NGC2623 








NGC2992 








NGC34 










1 1+ 1 - 4 
1 1± 0.8 


' " 




NGC4138 


04± ?:3 




NGC4151 




5 - 7± !:o 


4.3±« 4 


NGC4168 




NGC424 








NGC4258 




7.6±° : « 

1.4± f 
1.0 




NGC4303 




0.3±°;| 


NGC4395 






NGC4472 






22 - 3± ?60 


NGC4477 








NGC449 






NGC4507 




1.7±g-| 




NGC4565 






NGC4639 








NGC4725 




0.3±0;3 




NGC4945 






NGC4968 








NGC5033 


14± l'.2 






NGC5252 








NGC526A 








NGC5273 








NGC5506 
NGC5643 


1 4± 10 


2.o±f;i 




NGC591 






NGC6552 








NGC7172 








NGC7212 








NGC7314 








NGC7479 








NGC7582 




°' 5± 0.3 


0.7±g;2 


NGC7674 




UGC1214 




1 - 5± 0-9 




UGC2456 




12 ±o:8 


UGC2608 






UGC4203 








UGC6527 








UGC8621 








UM625 









Table 9. Best-fit fluxes for the RRCs discussed in this paper. Dots indicate lines not detected. 
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